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ABSTRACT 

We present a new grid of ionizing fluxes for O and Wolf-Rayet stars for use with evolu- 
tionary synthesis codes and single star Hn region analyses. A total of 230 expanding, 
non-LTE, line-blanketed model atmospheres have been calculated for five metallicities 
(0.05, 0.2, 0.4, 1 and 2Z Q ) using the WM-basic code of Pauldrach et al. (2001) for 
O stars and the CMFGEN code of Hillier & Miller (1998) for W-R stars. The stellar 
wind parameters are scaled with metallicity for both O and W-R stars. We compare 
the ionizing fluxes of the new models with the CoStar models of Schaerer & de Koter 
(1997) and the pure helium W-R models of Schmutz, Leitherer & Gruenwald (1992). 
We find significant differences, particularly above 54 eV, where the emergent flux is 
determined by the wind density as a function of metallicity. The new models have 
lower ionizing fluxes in the He I continuum with important implications for nebular 
line ratios. 

We incorporate the new models into the evolutionary synthesis code Starburst99 
(Leitherer et al. 1999) and compare the ionizing outputs for an instantaneous burst 
and continuous star formation with the work of Schaerer & Vacca (1998; SV98) and 
Leitherer et al. (1999). The changes in the output ionizing fluxes as a function of age 
are dramatic. We find that, in contrast to previous studies, nebular Hen A4686 will 
be at, or just below, the detection limit in low metallicity starbursts during the W-R 
phase. The new models have lower fluxes in the Hei continuum for Z ^ 0.4 Z Q and 
ages ^ 7 Myr because of the increased line blanketing. 

We test the accuracy of the new model atmosphere grid by constructing photoion- 
ization models for simple Hn regions, and assessing the impact of the new ionizing 
fluxes on important nebular diagnostic line ratios. For the case of an H n region where 
the ionizing flux is given by the WM-basic dwarf O star grid, we show that He I 
A5786/H/3 decreases between 1 and 2 Z in a similar manner to observations (e.g. 
Bresolin et al. 1999). We find that this decline is caused by the increased effect of line 
blanketing above solar metallicity. We therefore suggest that a lowering of the upper 
mass limit at high abundances is not required to explain the diminishing strength 
of He I A5786/H/?, as has been suggested in the past (e.g. Shields & Tinsley 1976; 
Bresolin et al. 1999). For an Hn region where the ionizing flux is provided by an 
instantaneous burst of total mass 10 6 M Q , we plot the softness parameter rf against 
the abundance indicator R23 for ages of 1-5 Myr. The new models are coincident with 
the observational data of Bresolin et al. (1999), particularly during the W-R phase, 
unlike the previous models of SV98 which generally over-predict the hardness of the 
ionizing radiation. 

The new model grid and upda ted Starburst99 code can be downloaded from 
http: / /www. star. ucl.ac.uk/starburst . 

Key words: stars: atmospheres - stars: mass loss - stars: Wolf-Rayet - H 11 regions 
- galaxies: starburst - galaxies: stellar content 
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observations at various wavelengths. In the satellite ultravio- 
let, the stellar wind spectral features of a population of mas- 
sive stars can be synthesized to provide information on the 
star formation rates, the slope of the initial mass function 
(IMF), and ages (e.g. Robert, Leitherer & Heckman 1993; 
Leitherer, Robert & Heckman 1995; de Mello, Leitherer & 
Heckman 2000). In the optical region, the properties of in- 
tegrated stellar populations are usually derived indirectly 
from their total radiative energy outputs using nebular diag- 
nostic line ratios (e.g. Garcia- Vargas, Bressan & Diaz 1995; 
Stasihska & Leitherer 1996; Stasihska, Schaerer & Leitherer 
2001). This method uses the theoretical ionizing fluxes from 
a population of massive stars as a function of age as input 
into a photoionization code. The ability of evolutionary syn- 
thesis models to predict correctly the properties of a young 
stellar population from nebular emission line ratios there- 
fore depends heavily on the accuracy of the evolutionary 
and atmospheric models developed for single massive stars. 

Early attempts to model stellar populations using evo- 
lutionary synthesis coupled with photoionization codes re- 
lied mainly on Kurucz (1992) plane-parallel LTE model at- 
mospheres (e.g. Garcia- Vargas et al. 1995). Gabler et al. 
(1989), however, showed that the presence of a stellar wind 
has a significant effect on the emergent ionizing flux in the 
neutral and ionized helium continua of O stars. Non-LTE 
effects depopulate the ground state of Hell, leading to a de- 
crease in the bound-free opacity above 54 eV, and hence a 
larger flux in the Hell continuum by up to 3-6 orders of 
magnitude. The emergent spectrum is flattened in the re- 
gion of the He I continuum and has a higher flux due to the 
presence of the stellar wind. 

Stasihska & Leitherer (1996) were the first to use ex- 
panding non-LTE atmospheres for stars with strong winds 
to construct photoionization models for evolving starbursts. 
They used the grid of pure helium, unblanketed non-LTE 
W-R atmospheres from Schmutz, Leitherer & Gruenwald 
(1992) to represent evolved stars with strong winds, and 
Kurucz (1992) models for hot stars close to the main se- 
quence. Line-blanketed, expanding non-LTE models cover- 
ing the main sequence evolution of O stars were first intro- 
duced by Stasihska & Schaerer (1997) who studied the effect 
of using the CoStar models of Schaerer & de Koter (1997) 
on the ionization of single star Hn regions. They found 
that higher ionic ratios were obtained in comparison to Ku- 
rucz models with the same stellar temperatures. Schaerer 
& Vacca (1998) used the CoStar and Schmutz et al. (1992) 
grids to construct evolutionary synthesis models for young 
starbursts. They predicted strong nebular Hell A4686 in low 
metallicity starbursts containing W-R stars. The same evo- 
lutionary synthesis models were used by Stasihska, Schaerer 
& Leitherer (2001) combined with photoionization models 
to analyse the emission line properties of H II galaxies. 

One major question arises: how realistic are the ioniz- 
ing fluxes being used in the evolutionary synthesis studies 
outlined above? To address this, there have been numer- 
ous studies aimed at empirically testing the accuracy of hot 
star model atmospheres by analysing H II regions containing 
single stars with well defined spectral types and effective 
temperatures. Schaerer (2000) reviews the success of non- 
LTE model atmospheres in reproducing nebular diagnostic 
line ratios. Esteban et al. (1993) examined the accuracy of 
the W-R grids by photoionization modelling of nebulae ion- 



ized by single W-R stars. They found that the lack of line 
blanketing was most important for the coolest W-R stars. 
In recent years, new computational techniques have allowed 
line-blanketed W-R atmospheres to be calculated for a few 
stars (e.g. Schmutz 1997; de Koter, Heap & Hubeny 1997; 
Hillier & Miller 1998). Crowther et al. (1999) tested these 
models by seeking to find a consistent model that reproduced 
the stellar and nebular parameters of a cool WN star. They 
found that line blanketing plays a significant role in mod- 
ifying the W-R ionizing output although the models of de 
Koter et al. (1997) and Hillier & Miller (1998) predict quite 
different ionizing flux distributions below the He I edge. 

Bresolin, Kennicutt & Garnett (1999) modelled extra- 
galactic Hn regions using Kurucz and non-LTE models. 
They find that the temperatures of the ionizing stars de- 
crease with increasing metallicity, and suggest that this can 
be explained by lowering the upper mass limit for star for- 
mation. In a further study, Kennicutt et al. (2000) have anal- 
ysed a large sample of H n regions containing stars of known 
spectral types and effective temperatures. They confirm the 
dependence of stellar temperature on metallicity, although 
they note that this relationship depends on the correctness 
of the input ionizing fluxes, and particularly the amount of 
line blanketing. Oey et al. (2000) have presented a detailed 
comparison of spatially resolved H n region spectra to pho- 
toionization models to test how well the CoStar models and 
Schmutz et al. (1992) W-R models reproduce the observed 
nebular line ratios. They find that overall the agreement is 
within 0.2 dex but the nebular models appear to be too hot 
by ~ 1000 K, and suggest that the ionizing flux distribution 
may be too hard in the 41-54 eV range. On the other hand, 
Kewley et al. (2001) model a large sample of infrared star- 
burst galaxies, and find that the Schmutz et al. (1992) W-R 
atmospheres do not produce sufficient flux between 13.6- 
54 eV to match their observations. 

It is clear that the use of expanding non-LTE model 
atmospheres in the analysis of single and unresolved H n re- 
gions is a major improvement over using static LTE models 
to represent hot massive stars. The empirical studies de- 
scribed above show that, in addition, it is essential that the 
atmospheres are line-blanketed, particularly for studies of 
young stellar populations at solar or higher metallicities. 
The CoStar models used in more recent studies do incorpo- 
rate line blanketing but the lack of any line blanketing in the 
W-R atmosphere grid of Schmutz et al. (1992) is a serious 
deficiency. With the recent advances in computing and the 
development of large codes to calculate expanding, non-LTE 
line blanketed atmospheres, it is now feasible to compute a 
grid of realistic ionizing fluxes for OB and W-R stars. In 
this paper, we present such a grid for five metallicities from 
0.05-2 Z for use with the Starburst99 (Leitherer et al. 1999) 
evolutionary synthesis code and analyses of single star Hn 
regions. In Section |j| we present the details of the model at- 
mosphere grid, and in Section |^, we compare the predicted 
ionizing fluxes with previous single star models used in syn- 
thesis codes. In Section ^, we incorporate the new grid into 
Starburst99, and in Section |j| we examine the changes in 
the ionizing fluxes as a function of age for an instantaneous 
burst and continuous star formation. In Section ^, we use 
the photoionization code CLOUDY (Ferland 2002) to assess 
the differences in the nebular diagnostic line ratios of Hn 
regions ionized by single O stars and synthetic clusters as 
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a function of age. We discuss our results in Section [7| and 
present the conclusions in Section H. 



2 THE MODEL ATMOSPHERE GRID 

We have sought to compute a grid of expanding, non-LTE, 
line-blanketed model atmospheres that cover the entire up- 
per H-R diagram for massive stars with stellar winds. To this 
end, we have used two model atmosphere codes: the WM- 
basic code of Pauldrach et al. (2001) for O and early B stars; 
and the cmfgen code of Hillier and Miller (1998) for W-R 
stars. In total, we have calculated 230 model atmospheres 
using a 350 MHz Pentium II PC for WM-basic and a 1.3 
GHz Pentium IV PC for cmfgen. Each O and W-R model 
calculation typically took 3 and 9 hours to complete respec- 
tively. We have calculated the models for five metallicities: 
0.05, 0.2, 0.4, 1, and 2Zq, as defined by the evolutionary 
tracks of Meynet et al. (1994). 

2.1 O star model atmospheres 

The current status of modelling O stars is reviewed by Ku- 
dritzki & Puis (2000). To calculate realistic ionizing fluxes, 
the effects of line blocking and blanketing have to be in- 
cluded in any expanding non-LTE model atmosphere code. 
This makes the calculation complex and computer-intensive 
because in an expanding atmosphere, the frequency range 
that can be blocked by a single spectral line is increased, and 
thus thousands of spectral lines from different ions overlap. 
Schaerer & de Koter (1997) used a Monte-Carlo approach 
(Schmutz 1991) which has the advantage that the influence 
of millions of lines on the emergent spectrum can be cal- 
culated, but the disadvantage that line broadening is not 
included, and thus the line blanketing in the low velocity 
part of the wind is underestimated. The WM-basic code of 
Pauldrach et al. (2001) overcomes these drawbacks by using 
an opacity sampling technique and a consistent treatment 
of line blocking and blanketing to produce synthetic high 
resolution spectra. We have therefore used the WM-basic 
code of Pauldrach et al. (2001) because of its fuller treat- 
ment of line blanketing. We discuss the differences in the 
predicted i oniz ing fluxes of the CoStar and WM-basic codes 
in Section 3.1. 

In Tables 1 and 2, we present the model grid we have 
computed for O and early B dwarfs and supergiants. We 
defined the fundamental parameters (effective temperature 
Tefs; gravity log g; and photospheric radius R t ) of the stellar 
atmospheres using the following approach. A revised abso- 
lute visual magnitude, My- spectral type scale for O dwarfs 
and supergiants was determined for Galactic and LMC stars, 
following the approach and references cited by Crowther & 
Dessart (1998) in their study of early O stars. A fit was made 
to individual absolute magnitudes as illustrated in Fig. |l|. 

Next, we adopted the revised effective temperature vs. 
spectral type calibration of Crowther (1998) for supergiants 
and dwarfs, which was updated relative to Vacca, Garmany 
& Shull (1996) to take account of more recent hydrostatic re- 
sults which neglect stellar winds and line blanketing. There 
is increasing evidence that the temperature scale for O 
dwarfs and especially supergiants based on hydrostatic mod- 
els is too high (Martins, Schaerer & Hillier 2002; Crowther 
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Figure 1. The adopted absolute magnitude My - spectral type 
calibrations for dwarf (solid), giant (dashed) and supergiant (dot- 
ted) O stars. The symbols indicate the mean values for each spec- 
tral type for dwarfs and supergiants. The typical error on these 
values is ±0.5 mag. 

et al. 2002b). Fortunately, the effect for O dwarfs is rela- 
tively small ~ 5 kK, but future modifications will certainly 
be necessary. 

We defined an effective temperature grid to span the 
domain of O and early B stars from 25 000-51 000 K, spaced 
at typically 0.03 dex intervals to give 11 models cover- 
ing the spectral classes from 03 to B1.5. Using the T e g- 
bolometric correction calibration of Vacca et al. (1996) we 
determined bolometric magnitudes and photospheric radii 
for each model. Since a discrepancy remains between masses 
derived from spectroscopic and evolutionary models (Her- 
rero et al. 1992), we adopted (logarithmic) mean surface 
gravities from Vacca et al. in all cases. 

The appropriate terminal velocity Voo for each model 
at solar metallicity was determined using the Voo -spectral 
type calibrations of Prinja, Barlow & Howarth (1990) and 
Lamers, Snow & Lindholm (1995). The mass loss rate M was 
derived from the wind momentum-luminosity relationship 
given by Kudritzki & Puis (2000) 

log [Mvoo (R,/R e ) - 5 ] = log Do + x log ^~ 

where log L> and x = 20.69,1.51 (OI); 19.87, 1.57 (OV); 
and 21.24, 1.34 (early BI). Since O star winds are driven 
through the absorption of photospheric photons by extreme 
ultraviolet (EUV) metal lines, mass loss rates and termi- 
nal velocities must depend on metallicity. Radiatively-driven 
wind theory (e.g. Kudritzki, Pauldrach & Puis 1987) pre- 
dicts that mass loss should scale with metallicity with an ex- 
ponent of 0.5-1.0. Recently Vink, de Koter & Lamers (2001) 
predict an exponent of 0.7 using Monte Carlo modelling. We 
have chosen to follow the prescription of Leitherer, Robert 
& Drissen (1992) in scaling both M and with metallicity 
by using power law exponents of 0.8 and 0.13 respectively 
i.e. 

Mz = Z°- s M ; Voo ,z = Z 013 Vo o, e 

The scaled values for M and v x are given for each metal- 



Table 1: OB Dwarf Model Atmosphere Grid 



Z 0.2 Z 0.4 Z 0.05 Z Q 2.0 Z q 



Model 




logs 


R* 


SpT 


log M 


Voo 


log M 


Voo 


log M 


Voo 


log M 


Voo 


log M 


Voo 


Hot. 


(kK) 


(cms ) 
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Qo 


Qi Q2 


Qo 


Qi Q2 


Qo 


Qi Q2 


Qo 


Qi Q2 


Qo 


Qi Q2 


OB#l 


50.0 


4.00 


9.8 


03V 


-5.85 


3150 


-6.41 


2550 


-6.17 


2790 


-6.89 


2130 


-5.61 


3440 












49.5 


48.8 45.8 


49.5 


49.0 45.6 


49.5 


48.9 45.7 


49.6 


49.1 45.4 


49.5 


48.8 45.6 


OB#2 


45.7 


4.00 


10.4 


04V 


-6.08 


2950 


-6.63 


2390 


-6.39 


2610 


-7.12 


1990 


-5.83 


3220 












49.4 


48.7 45.0 


49.4 


48.7 45.2 


49.4 


48.7 45.2 


49.4 


48.8 45.0 


49.4 


48.6 44.8 


OB#3 


42.6 


4.00 


10.5 


05V 


-6.24 


2870 


-6.80 


2330 


-6.55 


2550 


-7.28 


1940 


-6.00 


3140 












49.2 


48.5 44.4 


49.2 


48.6 44.2 


49.2 


48.5 44.3 


49.2 


48.6 44.2 


49.2 


48.5 44.5 


OB#4 


40.0 


4.00 


10.6 


07V 


-6.34 


2500 


-6.90 


2020 


-6.66 


2210 


-7.38 


1690 


-6.10 


2730 












49.0 


48.2 43.8 


49.0 


48.3 43.4 


49.0 


48.3 43.5 


49.0 


48.4 43.6 


49.0 


48.1 37.7 


OB#5 


37.2 


4.00 


10.5 


07.5V 


-6.47 


2100 


-7.03 


1700 


-6.79 


1860 


-7.51 


1420 


-6.23 


2290 












48.7 


47.6 37.3 


48.8 


47.8 42.3 


48.8 


47.8 42.4 


48.7 


47.8 42.5 


48.7 


47.4 38.5 


OB#6 


34.6 


4.00 


10.5 


08V 


-6.64 


1950 


-7.20 


1580 


-6.96 


1730 


-7.68 


1320 


-6.40 


2130 












48.5 


46.6 36.3 


48.4 


46.9 36.7 


48.4 


46.8 36.7 


48.4 


46.7 38.7 


48.4 


46.7 35.4 


OB#7 


32.3 


4.00 


10.1 


09V 


-6.74 


1500 


-7.30 


1210 


-7.06 


1330 


-7.79 


1010 


-6.50 


1640 












47.9 


45.7 34.3 


48.0 


45.5 34.8 


47.9 


45.5 34.9 


48.0 


45.6 34.8 


48.0 


45.6 33.6 


OB#8 


30.2 


4.00 


9.5 


B0V 


-6.92 


1200 


-7.48 


970 


-7.24 


1060 


-7.96 


810 


-6.68 


1310 












47.4 


44.8 32.7 


47.4 


44.6 33.3 


47.3 


44.8 33.4 


47.2 


44.5 33.2 


47.4 


44.9 32.7 


OB#9 


28.1 


4.00 


8.6 


B0.5V 


-6.74 


800 


-7.30 


640 


-7.06 


710 


-7.79 


540 


-6.50 


870 












47.0 


44.5 31.5 


46.8 


44.2 32.4 


46.8 


44.3 32.2 


46.8 


43.9 32.1 


46.9 


44.5 31.1 


OB#10 


26.3 


4.00 


8.1 


B1V 


-6.89 


700 


-7.45 


560 


-7.20 


620 


-7.93 


470 


-6.65 


760 












46.5 


43.9 30.5 


46.3 


43.7 31.6 


46.4 


43.8 31.4 


46.2 


43.3 30.9 


46.4 


44.1 


OB#ll 


25.0 


4.00 


7.1 


B1.5V 


-7.08 


600 


-7.63 


480 


-7.39 


530 


-8.12 


400 


-6.83 


650 












46.1 


43.4 


45.9 


43.2 30.8 


46.1 


43.3 30.9 


45.8 


42.9 


46.1 


43.7 



Notes: Mass loss rates M and terminal velocities Voo are in units of M yr 1 and kms 1 respectively. Qo, Qi and Q2 are on a log scale and represent the photon luminosity 
(s _1 ) in the Hi, He I and Hell ionizing continua. 



Table 2: OB Supergiant Model Atmosphere Grid 
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Q 
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Q 
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Qo 


Q 
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15.2 
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2790 




-5.63 




2130 




-4.34 
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50.0 


49 
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50.0 
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.5 46. 


.7 


50.0 


49 


.4 46 


.6 


50.0 


49. 


5 46 


.4 


50.0 


49 
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OB#13 


45.7 


3.73 


18.4 
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2320 


-5.28 




1880 




-5.04 




2060 
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1570 




-4.48 
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49.9 


49 


.0 38.2 


49.9 


49. 


.2 45. 


.6 


49.9 


49 


.2 38 


.9 


49.9 


49. 


3 45 


.9 


50.0 


48 


.9 37.2 


OB#14 


42.6 


3.67 


19.7 


051 


-4.92 




2300 


-5.48 




1860 




-5.24 




2040 




-5.96 




1550 




-4.68 




2510 












49.8 


48 


.9 38.3 


49.8 


49. 


.1 45 


A 


49.9 


49 


.1 45 


.2 


49.8 


49. 


2 45 


.4 


49.8 


48 


.7 37.6 


OB#15 


40.0 


3.51 


21.5 


061 


-5.01 




2180 


-5.57 




1760 




-5.33 




1930 




-6.05 




1470 




-4.77 




2380 












49.8 


48. 


.8 38.0 


49.8 


49. 


.0 44. 


.5 


49.8 


49 


.0 38 


.1 


49.8 


49. 


.1 38 


.2 


49.7 


48 


.6 36.9 


OB#16 


37.2 


3.40 


23.6 


07.51 


-5.07 




1980 


-5.62 




1600 




-5.38 




1750 




-6.11 




1340 




-4.82 




2160 












49.6 


48 


.5 37.2 


49.7 


48. 


.9 37. 


.6 


49.7 


48. 


.8 37 


.7 


49.7 


48. 
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.4 


49.6 


48 
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3.29 
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-6.23 
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45 


,7 


OB#21 


26.3 


2.99 


26.3 


BOI 


-6.08 




1400 


-6.64 




1130 




-6.40 




1240 




-7.12 




940 




-5.84 




1530 












47.8 


45 


.2 


47.5 


44. 


.8 31. 


.8 


47.6 


45 


.0 31 


.5 


47.6 


44. 


6 32 


.2 


47.8 


45 


.3 


OB#22 


25.0 


2.95 


26.9 


B0.5I 


-6.11 




1200 


-6.67 




970 




-6.43 




1060 




-7.15 




810 




-5.87 




1310 












47.6 


44 


.9 


47.1 


44 


.2 31. 


.5 


47.4 


44 


.6 




47.1 


44. 


2 




47.7 


45. 


.1 



Notes: Mass loss rates M and terminal velocities Voo are in units of M yr 1 and kms 1 respectively. Qo, Qi and Q2 are on a log scale and represent the photon luminosity 
(s _1 ) in the Hi, He I and Hell ionizing continua. 
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Table 3: WN Model Atmosphere Grid I. 

Z 0.2 Z Q 0.4 Z 0.05 Z 2.0 Z Q ™ 



Model 


3* 


R* 


Y 


log M 


Voo 


log M 




Voo 


log M 




Voo 


log M 


Voo 


log M 


Voo 


Ref. 


(kK) 


(Rq) 




Qo 


Qi Q2 


Qo 


Q 


1 Q2 


Qo 


Q 


1 Q2 


Qo 


Qi Q2 


Qo 


Qi Q2 


WN#1 


30.0 


20.3 


0.76 


-4.93 


1100 


-5.49 




890 


-5.25 




970 


-5.97 


740 


-4.69 


1200 










48.8 


45.0 29.5 


48.8 


46. 


.8 34.7 


48.8 


46. 


.6 34.4 


48.8 


46.8 35.3 


48.7 


44.3 28.0 


WN#2 


32.5 


17.2 


0.78 


-4.91 


1170 


-5.47 




940 


-5.23 




1030 


-5.95 


790 


-4.67 


1280 










49.0 


47.4 34.3 


49.0 


47. 


7 37.0 


49.0 


47. 


7 36.8 


49.0 


47.7 36.5 


48.8 


44.3 28.3 


WN#3 


35.0 


14.9 


0.80 


-4.89 


1240 


-5.45 




1000 


-5.21 




1100 


-5.93 


840 


-4.65 


1350 










49.1 


47.8 37.1 


49.1 


48. 


1 37.3 


49.1 


48. 


.1 37.6 


49.1 


48.2 36.9 


49.0 


47.1 34.5 


WN#4 


40.0 


11.4 


0.85 


-4.83 


1370 


-5.39 




1110 


-5.15 




1210 


-5.87 


920 


-4.59 


1490 










49.2 


48.1 38.1 


49.2 


48. 


4 37.7 


49.2 


48. 


3 37.7 


49.2 


48.5 37.7 


49.2 


47.7 34.3 


WN#5 


45.0 


9.0 


0.91 


-4.76 


1520 


-5.32 




1230 


-5.08 




1340 


-5.80 


1020 


-4.52 


1660 










49.3 


48.4 38.1 


49.3 


48. 


6 38.2 


49.3 


48. 


6 38.1 


49.3 


48.7 36.9 


49.2 


47.9 35.2 


WN#6 


50.0 


7.3 


0.94 


-4.73 


1670 


-5.29 




1350 


-5.05 




1480 


-5.77 


1130 


-4.49 


1820 










49.3 


48.6 38.2 


49.3 


48. 


.8 38.6 


49.3 


48. 


.8 38.6 


49.3 


48.8 38.5 


49.3 


48.3 38.5 


WN#7 


60.0 


5.1 


0.98 


-4.69 


1960 


-5.25 




1580 


-5.01 




1730 


-5.73 


1320 


-4.45 


2140 










49.4 


48.9 38.4 


49.4 


49. 


.0 39.1 


49.4 


49. 


.0 39.2 


49.4 


49.0 38.9 


49.3 


48.6 38.6 


WN#8 


70.0 


3.7 


0.98 


-4.69 


2270 


-5.25 




1840 


-5.01 




2010 


-5.73 


1530 


-4.45 


2480 










49.4 


49.0 39.6 


49.4 


49. 


.1 39.8 


49.4 


49. 


1 39.8 


49.3 


49.0 39.8 


49.4 


48.8 38.6 


WN#9 


80.0 


2.8 


0.98 


-4.69 


2550 


-5.25 




2060 


-5.01 




2260 


-5.73 


1720 


-4.45 


2790 










49.4 


49.1 40.1 


49.3 


49. 


.1 40.4 


49.4 


49 


.1 40.3 


49.3 


49.1 47.5 


49.4 


48.9 38.9 


WN#10 


90.0 


2.3 


0.98 


-4.69 


2850 


-5.25 




2310 


-5.01 




2520 


-5.73 


1930 


-4.45 


3110 










49.3 


49.1 40.6 


49.3 


49. 


.1 47.9 


49.3 


49 


.1 40.4 


49.3 


49.1 48.0 


49.3 


49.0 39.5 


WN#11 


100.0 


1.8 


0.98 


-4.69 


3120 


-5.25 




2530 


-5.01 




2760 


-5.73 


2110 


-4.45 


3410 










49.3 


49.1 41.2 


49.3 


49. 


.1 48.3 


49.3 


49 


.1 48.1 


49.3 


49.1 48.3 


49.3 


49.0 40.4 


WN#12 


120.0 


1.3 


0.98 


-4.69 


3700 


-5.25 




3000 


-5.01 




3300 


-5.73 


2500 


-4.45 


4050 










49.3 


49.1 49.2 


49.2 


49. 


.1 48.5 


49.3 


49 


.1 48.5 


49.2 


49.1 48.6 


49.3 


49.1 41.1 



Notes: Mass loss rates M and terminal velocities Voo are in units of M yr 1 and kms 1 respectively. Qo, Qi and Q2 are on a log scale and represent the photon luminosity 
(s _1 ) in the Hi, He I and Hen ionizing continua. 



Table 4: WC Model Atmosphere Grid 



0.2 Z Q 0.4 Z 0.05 Z 2.0 Z 



Model 


7* 


R* 


Y 


log M 




log M 




Voo 


log M 


Voo 


log M 


Voo 


log M 


Voo 


Ref. 


(kK) 


(R©) 




Qo 


Qi Q2 


Qo 


Q 


1 Q2 


Qo 


Qi Q2 


Qo 


Qi Q2 


Qo 


Qi Q2 


WC#1 


40.0 


9.3 


0.56 


-4.72 


1250 


-5.28 




1010 


-5.04 


1100 


-5.76 


840 


-4.48 


1360 










48.9 


41.8 26.1 


49.0 


47. 


.3 33.0 


49.0 


47.6 34.0 


48.9 


46.9 31.4 


48.7 


42.1 25.9 


WC#2 


43.0 


8.0 


0.56 


-4.72 


1380 


-5.28 




1110 


-5.04 


1220 


-5.76 


930 


-4.48 


1510 










49.0 


42.1 27.5 


49.1 


47. 


.7 34.8 


49.1 


47.9 35.4 


49.0 


47.5 33.4 


48.8 


39.6 23.8 


WC#3 


46.0 


7.0 


0.56 


-4.72 


1510 


-5.28 




1220 


-5.04 


1340 


-5.76 


1020 


-4.48 


1650 










49.0 


46.7 31.3 


49.1 


48. 


.0 36.1 


49.1 


47.9 34.8 


49.1 


48.1 36.1 


48.9 


42.1 26.6 


WC#4 


50.0 


5.9 


0.56 


-4.72 


1650 


-5.28 




1330 


-5.04 


1460 


-5.76 


1110 


-4.48 


1800 










49.1 


47.5 33.1 


49.1 


48. 


3 37.0 


49.1 


48.1 36.6 


49.1 


48.4 35.3 


49.0 


42.2 27.3 


WC#5 


55.0 


4.9 


0.56 


-4.72 


1870 


-5.28 




1510 


-5.04 


1660 


-5.76 


1260 


-4.48 


2040 










49.1 


47.9 31.3 


49.2 


48. 


5 37.3 


49.2 


48.4 37.4 


49.2 


48.6 40.3 


49.0 


42.4 28.5 


WC#6 


60.0 


4.1 


0.56 


-4.72 


2100 


-5.28 




1700 


-5.04 


1860 


-5.76 


1420 


-4.48 


2290 










49.1 


48.1 36.7 


49.2 


48. 


7 37.6 


49.2 


48.6 37.5 


49.2 


48.7 39.9 


49.0 


47.9 34.8 


WC#7 


70.0 


3.0 


0.56 


-4.72 


2450 


-5.28 




1980 


-5.04 


2170 


-5.76 


1650 


-4.48 


2680 










49.1 


48.5 37.9 


49.2 


48. 


.8 38.3 


49.2 


48.7 37.9 


49.2 


48.8 40.3 


49.1 


47.9 32.9 


WC#8 


80.0 


2.3 


0.56 


-4.72 


2800 


-5.28 




2270 


-5.04 


2480 


-5.76 


1890 


-4.48 


3060 










49.2 


48.6 38.2 


49.2 


48. 


.9 39.9 


49.2 


48.8 38.6 


49.2 


48.9 40.3 


49.1 


48.2 37.3 


WC#9 


90.0 


1.8 


0.56 


-4.72 


3200 


-5.28 




2590 


-5.04 


2840 


-5.76 


2160 


-4.48 


3500 










49.2 


48.7 38.4 


49.2 


48. 


.9 42.6 


49.2 


48.9 39.3 


49.2 


48.9 46.9 


49.1 


48.4 38.3 


WC#10 


100.0 


1.5 


0.56 


-4.72 


3600 


-5.28 




2920 


-5.04 


3190 


-5.76 


2430 


-4.48 


3930 










49.2 


48.8 38.8 


49.2 


48. 


.9 40.7 


49.2 


48.9 40.3 


49.1 


48.9 47.7 


49.1 


48.5 35.6 


WC#11 


120.0 


1.0 


0.56 


-4.72 


4400 


-5.28 




3570 


-5.04 


3900 


-5.76 


3000 


-4.48 


4810 










49.2 


48.8 39.4 


49.1 


48. 


9 48.1 


49.2 


48.9 47.6 


49.1 


48.9 48.2 


49.1 


48.6 38.8 


WC#12 


140.0 


0.8 


0.56 


-4.72 


5200 


-5.28 




4220 


-5.04 


4620 


-5.76 


3520 


-4.48 


5690 










49.1 


48.9 40.0 


49.1 


48. 


.9 48.8 


49.1 


48.9 48.2 


49.0 


48.9 48.4 


49.1 


48.7 39.1 



Notes: Mass loss rates M and terminal velocities Voo are in units of Mq yr 1 and kms 1 respectively. Qo, Qi and Q2 are on a log scale and represent the photon luminosity 
(s _1 ) in the Hi, He I and Hell ionizing continua. 
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licity in Tables 1 and 2. To compute the models, we used 
version 1.16 of the WM-basic code which includes atomic 
models of all the important ions of 149 ionization stages of 
26 elements from H to Zn. The wind structure was solved 
for consistently in the hydrodynamic part of WM-basic, cor- 
responding to a velocity law with exponent j3 ~ 0.8. Inter- 
vention was required in order that the required mass loss 
rate and terminal velocity were obtained, for the stellar in- 
put parameters, via "force multipliers" which describe the 
radiative line acceleration. We tuned the force multipliers 
for each model by iteration to achieve the wind parameters 
given in Tables 1 and 2 (M to within 5 per cent and v x 
to within ±100 kms -1 ). We have not included the effect of 
shocks in the stellar wind since shock characteristics are ill- 
defined across the O spectral range. The presence of shocks 
will increase the ionizing flux at EUV and X-ray wavelengths 
(see Pauldrach et al. 2001 for details). In Tables 1 and 2, we 
also provide the photon luminosity in the Hi (logQo), He I 
(logQi) and Hell (logQ2) ionizing continua for each model. 



2.2 W-R star model atmospheres 

The basic parameters of W-R stars have been substantially 
revised in recent years, due to the availability of increasingly 
sophisticated model atmospheres that are now able to repro- 
duce their complex emission line spectra (see e.g. reviews of 
Hillier (1999); Crowther (1999)). 

We have calculated a grid of expanding, non-LTE, line- 
blanketed WN and WC model atmospheres with the CMF- 
GEN code of Hillier & Miller (1998) using what we consider 
to be the most appropriate parameters for these stars. Line 
blanketing is accounted for in this code using the concept 
of super-levels, with reasonably sophisticated model atoms 
for the dominant ionization stages, namely H, He, C, N, O, 
Ne, Si, S, Fe in WN stars and He, C, O, Si, Fe in WC stars. 
Subsequent to the completion of this grid, computational 
improvements permitted many other elements and ioniza- 
tion stages to be included in CMFGEN, but tests indicate 
that these tend to be second order effects. 

In Tables 3 and 4, we present the parameters adopted in 
our WN and WC grids, parameterized by the stellar temper- 
ature T* , the stellar radius _R» (both defined at a Rosseland 
optical depth of ~ 10; see Section ^ for a discussion of W- 
R temperature definitions), and the helium mass fraction 
Y . To define these parameters, we used results from recent 
studies presenting line-blanketed model atmosphere analy- 
ses of W-R stars. For the WN class, our main references 
were Herald, Hillier & Schulte-Ladbeck (2001) and Hillier & 
Miller (1998), supplemented by Crowther & Smith (1997) 
for H/He ratios. For the WC class, we were guided by the 
analyses of Hillier & Miller (1999), Dessart et al. (2000) and 
Crowther et al. (2002a). Nugis & Lamers (2000) present a 
compilation of the stellar parameters for 64 Galactic W-R 
stars and provide a good discussion of the uncertainties. 

We defined a T* range of 30 000-120 000 K for the WNs 
and assigned a decreasing hydrogen content from 24 per cent 
(representing a cool late WN star) to 2 per cent (repre- 
senting a hot WNE star with T» > 60 000 K), plus CNO 
equilibrium values in all cases. For the WC stars, we chose 
a T« range of 40 000-140 000 K, and fixed C/He=0.2 and 
C/0=4 by number. We assumed representative luminosities 



of L = 3 x 10 5 L Q and 2 x 10 s L for the WN and WC models 
respectively. 

The winds of W-R stars are believed to be inhomoge- 
neous or clumped. Evidence for this comes from observa- 
tions of line profile variability (e.g. Moffat et al. 1988) and 
the over-prediction of the strength of the scattering wings 
in theoretical line profiles (Hillier 1991). The significance 
of clumping is that it leads to a downward revision in W- 
R mass loss rates although the exact amount is unclear, 
since line profile analyses give M/y/f where / is the vol- 
ume filling factor (e.g. Hillier & Miller 1999). In an analysis 
of the infrared and radio continuum fluxes for a large sam- 
ple of W-R stars, Nugis, Crowther & Willis (1998) derived 
dumping-corrected mass loss rates a factor of 3-5 lower than 
the smooth wind values. This study, however, used a radio 
sample which is biased towards higher mass loss rates, and 
thus the clumping factors should be viewed as upper lim- 
its (C. Leitherer, priv. comm.). Recent line profile analyses 
using CMFGEN find good consistency using clumping factors 
of / ~ 0.1 (e.g. Dessart et al. 2000; Crowther et al. 2002a), 
which correspond to mass loss rates that are a factor of ~ 3 
times lower than homogeneous models. 

To derive representative mass loss rates at solar metal- 
licity, we used the mass loss-luminosity-chemical composi- 
tion relationships given by Nugis & Lamers (2000) which are 
based on the dumping-corrected mass loss rates of Nugis et 
al. (1998). Specifically for the WN stars, we used the rela- 
tionship 

log M = -13.60 + 1.63 log + 2.22 log Y 

Lq 

which agrees well with results of line-blanketed analyses 
(Crowther 2001). For WC stars, Nugis & Lamers (2000) find 
that their mass loss rates do show an increase with luminos- 
ity but the scatter is large. Since we have chosen constant 
values of L and Y for all WC subtypes we used the WC 
linear regression relation of Nugis & Lamers (2000) to give 
a fixed mass loss rate of log M = —4.72 at solar metallicity. 
To assign terminal velocities, we determined approximate 
spectral types and used the Uoo-spectral type calibrations 
of Prinja et al. (1990). 

It is not known whether the wind parameters of W-R 
stars scale with metallicity because the driving mechanism 
of their winds has not been conclusively determined, and 
their atmospheres are chemically enriched in the products of 
the CNO-cycle (WNs) and He-burning (WCs). In the past, 
it has usually been assumed that there is no scaling with 
metallicity. 

Observationally, the parameters of Galactic and LMC 
WN stars do not show a clear metallicity dependence 
(Hamann & Koesterke 2000) although the factor of ~ two 
change in metallicity may not be sufficient to reveal sig- 
nificant variations within the observational and modelling 
uncertainties. Crowther (2000) analysed one SMC WN star 
and again found no differences relative to comparable stars 
in the LMC and Galaxy, although the spread in properties 
was large. More recently, Crowther et al. (2002a) have stud- 
ied a sample of LMC WC4 stars and derived wind densities 
which were ~0.2 dex lower than Galactic WC stars anal- 
ysed in the same manner (Dessart et al. 2000). Crowther et 
al. concluded that a metallicity dependence for W-R stars 
naturally explains the well known difference between the 
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Abundances M Scaling 

Zq Zq 

0.05 Z Q 0.05 Z 

0.05 Zq Z 



2.2 2.4 2.6 2.8 3.0 

log A (A) 

Figure 2. The flux distribution of model WC#10 (T„ = 
100 000 K) for three different assumptions about the abundances 
and mass loss scaling: (a) abundances and mass loss rate scaled 
to solar (black, solid); (b) abundances and mass loss rate scaled 
to 0.05 Zq (red, dashed); and (c) 0.05 Zq abundances but mass 
loss rate scaled to solar (blue, dotted). Only model (b) with wind 
parameters scaled to 0.05 Zq has a significant flux below the He+ 
edge at 228 A because of the reduced wind density. The flux is in 
units of ergs cm -2 s — 1 A -1 measured at the stellar surface. 



Galactic and LMC WC populations, since the classification 
diagnostic C ill A5696 is extremely mass-loss dependent. 

Recent theoretical work also indicates that the strength 
of W-R winds should be metallicity dependent. Although the 
precise wind driving mechanism is still unsolved, it is gener- 
ally believed that radiation pressure with multiple scatter- 
ing is responsible. The reduction in mass loss rates through 
clumping has alleviated the momentum problem but the mo- 
mentum transfer efficiencies still lie in the range 2-20 (Nugis 
& Lamers 2000). Heavy metals are critical for initiating the 
outflow, even for WC stars, since the high ions of C and 
O do not possess enough driving lines in the correct energy 
range (Crowther et al. 2002a). Therefore we have adopted 
the same power law expo nent scaling as used previously for 
the O star models (Sect. 2.1). The weak metallicity effect 



identified by Crowther et al. (2002a) is consistent with the 
scaling in our approach. 

Observations show that nebular Hell A4686 emission 
is spatially associated with hot W-R stars in low metallic- 
ity environments (Garnett et al. 1991; Schaerer, Contini & 
Pinado 1999). Schmutz et al. (1992; hereafter SLG92) dis- 
cuss the fact that the energy output above 54.4 eV depends 
on the mass loss rate. In a W-R star with a strong wind, 
He ++ recombines at some point, whereas there is no recom- 
bination in a weak wind, and flux is emitted above 54.4 eV. 
In Fig. |^, we show the flux distribution of a WC model 
with T* = 100 000 K for three different assumptions about 
the line blanketing and wind density. It is apparent that for 
a solar (i.e. unsealed) mass loss rate, no ionizing photons 
escape below 228 A at Z = 0.05 and Z =7iq. Indeed, the 
emergent spectra are very similar, demonstrating that the 
reduced line blanketing at Z = 0.05 Zq has a small effect. 
In contrast, the model which has wind parameters scaled 
to Z — 0.05 Zq has a significant ionizing flux below the 
He + edge. We thus find that the wind density controls the 
transparency of the wind below 228 A at a given stellar tem- 
perature and the effect of line blanketing is negligible. The 



photon luminosities of the models are given in Tables 3 and 4 
in the form of log Qo, log Qi and log Q^. It should be noted 
that no models with 2 Zq abundances have significant Q-2 
fluxes. 

The mass loss rates in Tables 3 and 4 are corrected for 
clumping since they are derived from the Nugis & Lamers 
(2000) relationship. For consistency, we have introduced 
clumping into the CMFGEN model atmosphere calculations 
by adopting a volume filling factor / of 0.1 (Hillier & Miller 
1999). This has no effect on the emergent flux but does make 
a difference to the shape of the W-R line profiles. It is pos- 
sible to recover the unclumped mass loss rates by multi- 
plying the values in Tables 3 and 4 by l// 1 ^ 2 , i.e. VTO. 
For these calculations, we assume a standard velocity law of 
v(r) — i>oo(l — R/r) 13 with (3 = 1, and a Doppler line width 
of 200 kms -1 . 

Since the wind density is such a critical parameter, 
our use of the Nugis & Lamers (2000) mass-loss-luminosity 
calibration deserves special mention. As discussed above, 
this reproduces the average mass-loss rates for strong-lined 
Galactic W-R stars rather well. The weak-lined W-R stars, 
however, fall up to ~0.5 dex below this calibration. Some 
pathological stars have even weaker winds and are not ac- 
counted for in our calculations. Examples are HD 104994 
(WN3) which Crowther, Smith & Hillier (1995) predicted 
to emit strongly in the He + continuum, whilst nebular He II 
A4686 is observed in the Hn region G 2.4+1.4 which is as- 
sociated with Sand 4 (WO, Esteban et al. 1992). LMC W-R 
stars known to emit strongly below A228A also exist, such 
as Brey 2 (WN2, Garnett & Chu 1994). Fortunately, such 
W-R stars are relatively few in number, so it is probable 
that neglecting their contribution is acceptable for the large 
massive star population in a starburst. 



3 MODEL ATMOSPHERE COMPARISONS 
3.1 O star models 

The first expanding, line-blanketed, non-LTE grid of O star 
model atmospheres was introduced by Schaerer & de Koter 
(1997; hereafter SK97). These combined stellar structure 
and atmosphere or 'CoStar' models are used extensively in 
single star Hn region analyses (e.g. Oey et al. 2000) and 
spectral synthesis studies (e.g. Stasihska et al. 2001). The 
CoStar grid consists of 27 models covering the main se- 
quence evolution of 03 to early B spectral types at metal- 
licities of solar and 0.2 Zq. The mass loss rates are taken 
from Meynet et al. (1994) and are scaled with metallicity 
using an exponent of 0.5. To compare the WM-basic with 
the CoStar models, we have chosen a hot (T e ff = 43 kK; 
OB#14) early O supergiant model with a strong wind, and 
a cool (T e ff = 33 kK; OB#7) late O dwarf model with a weak 
wind. In Figs. ^ and ^ we show the comparisons at solar and 
0.2 Zq and also plot Kurucz (1992) LTE plane parallel mod- 
els for reference. The details of the model parameters are 
given in the figure captions. 

There are clearly significant differences in the emergent 
fluxes from the WM-basic and CoStar models. Considering 
first the O supergiant case, the WM-basic model has no flux 
below the He + edge at Zq in contrast to the CoStar model 
whereas at 0.2 Zq, both models have flux below 228 A. We 
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Figure 3. Comparison of the emergent fluxes for an early O su- 
pergiant from a WM-basic (black, solid), a CoStar (red, dashed) 
and a Kurucz model (blue, dotted) at solar metallicity and 0.2 Zq 
The model details for Zq are: WM-basic OB#14 T e g = 42.6 kK 
logg = 3.67; H* = 19.7R ; logM = -4.92; v x = 2300 kms" 1 
CoStar #E3 T eff = 42.6 kK; logg = 3.71; R* = 20.0 R© 
logM = -4.88; = 2538 kms" 1 ; and Kurucz T cS = 42.5kK 
logg = 5.0. The model parameters for 0.2 Zq are scaled as dis- 
cussed in the text. 



find again that the wind density is the critical parameter in 
determining the transparency below 228 A for the WM-basic 
models; a model with an un-scaled (i.e. solar) mass loss rate 
and 0.2 Zq abu ndances has no flux below this wavelength 
(cf. Section 2.2). The reason for the difference between the 
CoStar and WM-basic models at Zq is less clear because 
both models have the same wind density. Crowther et al. 
(1999) explore the differences in the ionizing fluxes predicted 
by CMFGEN and the ISA-wind code of de Koter, Heap & 
Hubeny (1997) which is used for the CoStar atmosphere 
calculations. They suggest that the most likely reason for 
the difference is that ISA-wind overestimates the far-UV flux 
because photon line-absorption and subsequent re-emission 
at longer wavelengths is not taken into account. 

For the late O dwarf case, we find that the WM-basic 
model has a lower flux in the He continuum compared to 
both the CoStar and Kurucz models. It is closer to the Ku- 
rucz model, as expected for cool dwarf stars with low mass 
loss rates. The lower flux compared to the Kurucz model 
is probably due to non-LTE effects producing deeper line 
cores in the blocking lines (Sellmaier et al. 1996). SK97 in- 
vestigated the excess flux in their dwarf models compared 
to plane parallel non-LTE models. They find that the tem- 
perature is too high in the region where the He continuum 
is formed and caution against the use of these models. 
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Figure 4. Comparison of the emergent fluxes for a late O dwarf 
from a WM-basic (black, solid), a CoStar (red, dashed), and a 
Kurucz model (blue, dotted) at solar metallicity and 0.2 Zq. The 
model details for Zq are: WM-basic OB#7 T eff = 32.3 kK; log g = 
4.00; R, = IO.IRq; logM = -6.74; v x = 1500 kms" 1 ; CoStar 
#A2 T cff = 33.3 kK; logg = 4.01; R, = 7.3 Rq; log M = -6.86; 
Voa = 2544 kms" 1 ; and Kurucz T cff = 32.0 kK; logg = 4.0. The 
model parameters for 0.2 Z are scaled as discussed in the text. 



3.2 W-R models 

As discussed in Section the W-R grid is based on current 
best estimates of wind parameters, temperatures, luminosi- 
ties and chemical compositions. These parameters differ sub- 
stantially from the pure helium W-R model grid of Schmutz 
et al. (1992) which is implemented in the Starburst99 evolu- 
tionary synthesis code. The SLG92 grid is defined by three 
parameters: the 'transformed radius' Rt (a measure of the 
inverse wind density) given by 

/ io-4 \ 2 / 3 

R-R 10 V °° 
1 * \dM/dt 2500 

the stellar or core effective temperature T t , and the veloc- 
ity law exponent (ft = 1 or 2). In Starburst99, SLG92 W-R 
atmospheres are connected to the evolutionary tracks by in- 
terpolating T» and Rt to match the hydrostatic effective 
temperature Th y d tabulated in Meynet et al. (1994) and the 
M and values of Leitherer et al. (1992) as a function 
of W-R type. As discussed in the next section, the interfac- 
ing of W-R atmospheres to evolutionary tracks is difficult 
because of the optically thick nature of the W-R winds. We 
have chosen to implement the CMFGEN W-R atmospheres by 
matching T* to 0.6Th y d + 0.4T 2 /3 where T 2 /3 represents the 
corrected temperature tabulated by Meynet et al. (1994). To 
make a meaningful comparison between the SLG92 and the 
new W-R grids, we compare models that would be selected 
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Figure 5. The emergent fluxes for CMFGEN WN models with 
T* = 45,60 and 90 kK for Z =Z Q (black, solid) and 0.2 Z Q (blue, 
dotted). For comparison, the SLG92 models are also plotted (red, 
dashed) that would be chosen in Starburst99 for the same evo- 
lutionary point at 4Myr in the HRD but corresponding to the 
hydrostatic core Meynet et al. temperatures of 56.2, 82.2 and 
130.1 kK. The transformed radii Rt of these three models are 
17.9, 8.4 and 1.7 R Q . 



in Starburst99 and our updated version to represent a W-R 
star at a given evolutionary point in the H-R diagram. 

In Fig. U we show three CMFGEN WN models with 
T* = 45 000, 60 000 and 90 000 K for metallicities of solar 
and 0.2 Zq. For comparison, we plot the equivalent SLG92 
models that would be chosen in Starburst99 at an age of 
4 Myr corresponding to T hyd = 56 200, 82 200 and 130 100 K. 
Significant differences are seen because of the different wind 
densities, our inclusion of line blanketing, and the lower tem- 
peratures we adopt to match the W-R atmospheres to the 
evolutionary tracks. For the 45 000 K model, the emergent 
flux below 504 A is much lower for the CMFGEN model at Zq 
compared to SLG92 because of line blanketing, as shown by 
the increased flux in the 0.2 Z model. The 60 000 K CMFGEN 
model has no flux below 228 A in contrast to the higher tem- 
perature SLG92 model, and the flux below 504 A is again re- 
duced. The large difference in the emergent flux below 228 A 
seen in Fig. || at 90 000 K at Zq between the CMFGEN and 
SLG92 mod el is due to the wind density effect discussed in 
Section \2:l\ since the 0.2 Zq CMFGEN model has emergent 
flux below 228 A. 

In Fig. we plot two CMFGEN WC models with tem- 
peratures of T„ = 60 000 and 100 000 K and compare them 
to the models that would be chosen in Starburst99 at an 
age of 5 Myr corresponding to T hyd = 78 200 and 126 100 K. 
The differences seen for the cooler case are similar to the 
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Figure 6. The emergent fluxes for CMFGEN WC models with 
T» = 60 and 100 kK for Z =Z (black, solid) and 0.2 Zq (blue, 
dotted). For comparison, the SLG92 models are also plotted (red, 
dashed) that would be chosen in Starburst99 for the same evo- 
lutionary point at 5 Myr in the HRD but corresponding to the 
hydrostatic Meynet et al. temperatures of 78.2 and 126.1 kK. The 
transformed radii Rt of these two models are 2.7 and 1.6 Rq. 

45 000 K WN model comparison in that the line blanketing 
at Zq produces a significant reduction in the flux in the 
He continuum. The differences seen at 100 000 K can be 
explained by the wind density determining the flux below 
228 A. 



4 INTEGRATION INTO STARBURST99 

Leitherer et al. (1999) present the evolutionary synthesis 
code Starburst99 which predicts the observable properties 
of galaxies undergoing active star formation, and is tailored 
to the analysis of massive star populations. The code is an 
improved version of that previously published by Leitherer 
& Heckman (1995). It is based on the evolutionary tracks of 
Meynet et al. (1994) for enhanced mass loss and the model 
atmosphere grid compiled by Lejeune et al. (1997), supple- 
mented by the pure helium W-R atmospheres of Schmutz et 
al. (1992). 

To integrate the new O star grid, we have simply re- 
placed the Lejeune et al. (1997) LTE models by the WM- 
basic models for effective temperatures above 25 000 K, re- 
mapped to 1221 points over the wavelength range log A = 
1.96-6.20 A, as required by the code. To ensure a smooth 
transition to the LTE models, we adopt the method of SV98: 
the WM-basic models are restricted to the T c g-log g domain 
defined by logp ^ 2.2 and log 5 < 5.71 x logT cff - 21.95. 
The ionizing fluxes of giants were determined by averaging 
the supergiant and dwarf models, and renormalising to con- 
serve the flux. Tests using real giant models show that this 
method is sufficiently accurate. We adopt the same method 
used by Leitherer et al. (1999) in choosing the best model 
atmosphere to represent an evolutionary point viz. selection 
using the nearest fit in log g and T a g and then interpolation, 
rescaling for T c s and R. 

The new W-R grid is based on current best estimates 
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Figure 7. Histogram plots of the temperature distribution for 
the WN and WC phases summed over the lifetime of the W- 
R phase as given by Starburst99 with the Meynet et al. (1994) 
evolutionary tracks at solar metallicity for an instantaneous burst 
with M = 10 6 Mq. For each W-R type, we show the temperature 
distributions corresponding to the hydrostatic temperature Th y( j 
(light grey), the corrected hydrostatic temperature T 2 /3 (dark 
grey) of Meynet et al. (1994), and the adopted T„ = 0.6T hyd + 
0.4T 2 /3 (black). The temperature bin sizes are 10 kK. 



of W-R temperatures and these are considerably lower than 
those of SLG92 and the evolutionary model temperatures. 
The problem of interfacing W-R atmospheres with evolu- 
tionary tracks is an unresolved issue which SLG92 discuss 
in some detail. It exists because the temperature given by 
evolutionary models Th y d refers to the hydrostatic core. 
Since W-R winds are optically thick, the observed radiation 
emerges at larger radii, and the temperature derived from 
observational analyses cannot be related to Th y d without in- 
ward extrapolation. To overcome this, Maeder (1990) and 
Meynet et al. (1994) correct the hydrostatic Thyd value by 
assuming a velocity law to give the radius corresponding to 
an optical depth of 2/3 and hence a lower temperature T2/3. 
This does not, however, solve the problem because the corre- 
sponding T2/3 temperatures derived from W-R model atmo- 
spheres analyses are generally close to 30 000 K irrespective 
of W-R subtype. As discussed by SLG92, W-R model at- 
mospheres are best characterized by the core temperature 
T* which refers to the radius where the optical depth is ~ 
10 (corresponding to a thermalization depth of unity for the 
continuum opacity which is dominated by electron scatter- 
ing)- 

The problem is therefore how to assign a W-R model 
atmosphere characterized by T, to an evolutionary model 
denned by two temperatures: Th y d and T2/3, with the former 



being too hot and the latter too cool. In Starburst99, an 
SLG92 W-R atmosphere is assigned by interpolating T* and 
Rt to match Th y d of the Meynet et al. (1994) evolutionary 
models. Since the Thyd values of Meynet et al. are higher 
than the new W-R grid values, we have assigned the W-R 
models to the tracks by matching T» to 0.6T hy d + 0.4T 2 / 3 . 
We arrived at this formalism by studying the distribution 
of temperatures in the W-R phase as given by Starburst99 
for solar metallicity. We simply binned the number of WN 
and WC stars into 10 kK temperature bins and summed over 
the lifetime of the W-R phase. The resulting WN and WC 
temperature distributions are shown in Fig. m 

For the WN stars, the temperature distribution for Th y d 
shows one broad peak at ~ 40 000 K and a second peak near 
120 000 K. The distribution for the WC stars is flatter and 
peaks at ~ 150 000 K. We first tried a simple average of the 
two Meynet et al. temperatures but this shifted the peak 
of the WN temperature distribution to a cool ~ 30 000 K. 
We then adjusted the bias between the two temperatures by 
trial and error and found that 0.6Th y d +0.4T 2 /3 gave a single 
temperature distribution for the WNs peaking at 40 000 K 
and extending to 100 000 K. The corresponding WC temper- 
ature distribution, as shown in Fig. m, peaks at 110 000 K. 
Ideally, the next step would be to compare the adopted T* 
distribution with observed W-R subtype vs. temperature 
distributions. Unfortunately there is insufficient data to do 
such a comparison. Nevertheless, our adopted T* values for 
WN and WC stars span the range occupied by Galactic and 
LMC W-R stars, as derived by recent analyses (e.g. Hillier 
& Miller 1998; Herald et al. 2001; Crowther et al. 2002a). 
We also note that the WN and WC T» peaks correspond to 
late WN and early WC stars. Studies of W-R galaxies (e.g. 
Guseva, Izotov & Thuan 2000) show that the composite W- 
R feature at ~ 4680 A is usually dominated by these types 
of W-R stars. 

Finally, we have adopted the same switching technique 
employed by Leitherer et al. (1999): W-R atmospheres are 
used when the surface temperature of the Meynet et al. 
(1994) tracks exceeds 25 000 K and the surface hydrogen 
content is less than 0.4. In addition, a WN or WC model 
is chosen depending on whether C/N< 10 (WN) or > 10 
(WC). 



5 EVOLUTIONARY SYNTHESIS 
COMPARISONS 

In Section ^, we compared the new O and W-R model at- 
mospheres with the grids of Schaerer & de Koter (1997) 
and Schmutz et al. (1992). We find that there are signifi- 
cant differences below 228 A at Zq (and higher) for the hot 
O supergiant and W-R models. We can thus expect large 
differences in the output ionizing fluxes of young massive 
star-forming regions, especially during the W-R phase. We 
compare the emergent flux distributions obtained with the 
new models integrated into Starburst99 to: 

(i) the standard version of Starburst99 (which we will re- 
fer to as SB99) which uses the W-R atmospheres of Schmutz 
et al. (1992) for stars with strong mass loss and the LTE 
models from the compilation of Lejeune et al. (1997) for all 
other stars. 
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Figure 8. Spectral energy distributions of the new models (left-hand side) compared with Starburst99 (right-hand side) at time intervals 
of 1, 3, 4, 5 and lOMyr for an instantaneous burst at Zq. 



(ii) a modified version of (i) with the CoStar models of 
Schaerer & de Koter (1997) replacing the LTE models for 
the main sequence evolution of massive stars. We have im- 
plemented these models as described in Schaerer & Vacca 
(1998) and use the solar metallicity grid for solar and above, 
and the 0.2 Zq grid for Z ^ 0.4 Zq. We refer to these models 
as SV98 and have checked that the output is the same as 
the actual SV98 models for the hot O star phase. There are 
small differences in the W-R phase resulting from the dif- 
ferent methods of interpolation used in the SB99 and SV98 
codes. 

We consider two types of star formation: an instantaneous 
burst and continuous star formation. We adopt a standard 
initial mass function (IMF) with a single Salpeter power 
law slope of 2.35 and lower and upper mass cut-offs of 1 and 
100 Mq . We use the enhanced mass loss tracks of Meynet et 
al. (1994) and calculate models for five metallicities: 0.05, 
0.2, 0.4, 1 and 2Zq. We do not include the contribution of 
the nebular continuum. 



5.1 Instantaneous burst 

We adopt a standard instantaneous burst model with a total 
mass of 10 6 M . In Fig. ^, we show the spectral energy dis- 
tributions obtained with the new model atmospheres com- 
pared with SB99 at time intervals of 1, 3, 4, 5 and lOMyr 
for Zq . The most dramatic differences occur below the He + 
edge at 228 A during the W-R phase at 3-5 Myr. The new 
line-blanketed models have negligible flux below 228 A at Zq 
in contrast to the very hard ionizing fluxes of the SLG92 at- 
mospheres. The detailed differences between the new models 
and previous work, and their dependence on metallicity are 
best seen by considering the ionizing fluxes and their ratios. 
In Fig. ^, we show the evolution of the photon luminos- 
ity from 1-8 Myr in the ionizing continua of Hi (logQo), 
Hel (logQi) and Hell (logQ2) for 0.2 and Zq. The evolu- 
tion of the hardness of the ionizing spectra log(Qi/Qo) and 



log(Q2/Qo) for all five metallicities is shown in Fig. [jlj. The 
variation with metallicity is primarily due to the changing 
W-R/O ratio. The relative number of W-R stars formed and 
the duration of the W-R phase increases with metallicity 
because the lower mass limit for their formation decreases 
(Meynet 1995). The W-R/O ratio changes from 0.025 to 0.36 
between Z = 0.05-2 Z Q at 3 Myr (SV98). 

The ionizing flux in the Lyman continuum Qo is de- 
termined mainly by hot main sequence stars. Fig. ^ shows 
that Qo is very similar to SV98 although the SB99 values 
are slightly lower because the Kurucz LTE models for hot 
O stars have smaller radii than the non-LTE O star models. 
The new models have lower ionizing fluxes below 504 A at 
ages of less than ~ 7 Myr for metallicities ^ 0.4 Zq. This 
is because W-R stars contribute an important fraction of 
the Hel ionizing flux and our new W-R models are line- 
blanketed compared to the SLG92 grid used for the SB99 
and SV98 models. The new O star models also have a lower 



ionizing flux in the Hel continuum (Sect. 3.1). The de- 



pendence of log(Qi/Qo) on metallicity is demonstrated in 
Fig. [h]. At 5 Myr, the ratio Qi/Qo is softer by a factor of 2 at 
2 Zq compared to SB99 and SV98 because of the effect of the 
W-R line-blanketed atmospheres, whereas at Z = 0.05 Zq, 
the models agree with SV98 because at this low metallicity, 
the main contributors to Qi are O supergiants because of 
the low W-R/O ratio. In general, at sub-solar metallicities 
for Qi/Qo, the new models agree more closely with SV98 
than SB99. 

The most pronounced revisions arc seen in Q2 during 
the W-R phase at 3-6 Myr in Fig. ^ for Z . The number of 
Hen ionizing photons emitted is now negligible compared 
to SV98 and SB99 because of the W-R wind density ef- 
fect (Sect. 2.2) and the conspicuous bump in the SV98 and 



SB99 models disappears. At metallicities below solar, Q2/Q0 
(Fig. |l"o| ) is softer in the W-R phase compared to SV98 and 
SB99 by, for example, a factor of ~ 25 at 0.4 Z Q . The W-R 
subtypes can be identified since we have used different mod- 
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Figure 9. The evolution of the photon luminosity for ages of 1—8 Myr in the ionizing continua of hydrogen (logQo)i He I (logQi) and 
Hen (logQ2) at 0.2 and Zq for an instantaneous burst at time steps of 0.1 Myr. The new models (black, solid) are compared to the 
ionizing fluxes of: Starburst99 (SB99, red, dotted) using the stellar library of Lejeune et al. (1997) and the Schmutz et al. (1992) W-R 
models; the Schaerer & Vacca (1998) models (SV98, blue, dashed) using CoStar and Schmutz et al. (1992) W-R atmospheres. 
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Figure 10. The evolution of the hardness of the ionizing spectra shown by the ratios of log(Qi/Qo) (lefthand side) and log(Q2/Qo) 
(righthand side) at all five metallicities for an instantaneous burst and ages of 1—8 Myr. The details of the models are given in Fig. bl 
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els for the WN and WC phases. At 0.4 Z , the first W-R 
stars to form at 3.5 Myr are the late WN stars, the descen- 
dants of the most massive O stars. These stars then evolve 
to WC stars, and finally, near the end of the W-R phase, 
the O stars near the minimum mass limit for W-R forma- 
tion, become WN stars. At Z = 0.05 Z Q , the W-R phase 
is short-lived because only the most massive O stars evolve 
into W-R stars. At all metallicities, the new O star models 
are softer compared to SV98 because the far-UV flux in the 
CoStar atm osph eres is probably overestimated (see discus- 
sion in Sect. |3.l|) . 

Overall, at metallicities of solar and higher, we find that 
the He n ionizing flux is negligible because of the high wind 
densities and significant line-blanketing. Hell ionizing pho- 
tons are produced during the W-R phase at metallicities 
below solar but Q2/Q0 is much softer by a factor of ~ 20 
compared to the SLG92 W-R models. In Section [?], we dis- 
cuss the presence of nebular Hell A4686 as an indicator of 
a significant W-R population and its strength as a function 
of metallicity. 



5.2 Continuous star formation 

For continuous star formation, we assume a star formation 
rate of 1 Mq yr _1 and the same IMF parameters adopted for 
the instantaneous case. In Fig. jn], we show the spectral en- 
ergy distributions for continuous star formation. In Figs, [l^ 
and[l^, the log Q values and their log ratios are shown. After 
~ 5 Myr, a balance is achieved between massive star birth 
and death and the ionizing fluxes become time independent. 

The differences shown between the various models in 
Fig. emphasize the differences between non-LTE and LTE 
models and blanketed vs. un-blanketed non-LTE models. 
Again, the predicted values of Qi are lower than SV98 be- 
cause of our use of line-blanketed W-R models and the lower 
ionizing fluxes of the WM-basic models. The difference in 
Q 2 for ages < 4 Myr between the SB99 and SV98 models 



is simply due to LTE vs. non-LTE models being used for 
O stars. The increase in Q2 when the first W-R stars are 
formed is much less marked compared to SB99 and SV98 
because of our use of line-blanketed W-R models. The low 
value of logQ2/Qo ~ —4 implies that nebular Hell will not 
be detectable at any metallicity from an integrated stellar 
population undergoing continuous star formation (Sect. Pzl) - 



6 IMPACT ON NEBULAR DIAGNOSTICS 

We now assess the impact of the new grid of model atmo- 
spheres on the ionization structure of H n regions by consid- 
ering single O star and synthetic cluster models as a function 
of age. To perform these tests, we have used the photoion- 
ization code CLOUDY (version 96, Ferland 2002). 

6.1 Single star Hll regions 

Many studies of the excitation of H n regions have suggested 
that the ionizing continuum becomes softer with increasing 
metal abundance (e.g. Shields 1974; Stasihska 1980; Vflchez 
& Pagel 1988). Shields & Tinsley (1976) first suggested 
that this relationship could be explained by a metallicity- 
dependent upper mass limit for star formation. Recently, 
Bresolin et al. (1999) and Kennicutt et al. (2000) have con- 
firmed the decrease in effective temperature with increasing 
metal abundance. They discuss the possibility of a lower 
upper mass limit operating at abundances above solar, al- 
though they caution that insufficient line blanketing in the 
model atmospheres at high abundances may be the root 
cause. We are now able to test this since our new O star 
grid contains fully line- blanketed models up to 2Zq. 

The nebular line ratio He I A5786/H/3 is a sensitive di- 
agnostic of the stellar temperature over the narrow range 
where helium is partially ionized, corresponding to w 35 000- 
40 000 K, but is insensitive to changes in the ionization pa- 
rameter U (Bresolin et al. 1999). Using the appropriate O 
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Figure 12. The evolution of the photon luminosity for ages of 1-10 Myr in the ionizing continua of hydrogen (logQo)i He I (logQi) 
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Figure 13. The evolution of the hardness of the ionizing spectra shown by the ratios of log(Qi/Qo) (lefthand side) and log(Q2/Qp 
(righthand side) at all five metallicities for continuous star formation and ages of 1—10 Myr. The details of the models are given in Fig. E 
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Figure 14. The predicted nebular He I A5786/H/3 ratio plotted 
against metallicity for effective temperatures between 37-45kK 
for WM-basic, CoStar and Kurucz O V model atmospheres. 



dwarf models listed in Table 1, we have calculated sim- 
ple dust-free, plane-parallel, ionization-bounded photoion- 
ization models with log U = —3.0 and a constant gas den- 
sity n = 50 cm -3 with a filling factor of unity. We adopt 
the abundances and depletion factors given by Dopita et al. 
(2000) for solar metallicity and scale these for other metal- 
licities. For the helium abundance as a function of metal- 
licity, we use the formula given by Dopita et al. (2000). To 
compare our results with the CoStar models, we have used 
the model grid available within CLOUDY corresponding to 
main sequence stars for metallicities of 0.2 Zq and Zq, and 
interpolated to match the temperatures of the WM-basic 
models. For Kurucz LTE models, we have used models that 
are the closest match to the WM-basic models in effective 
temperature and log g. 

In Fig. [jj, we plot He I A5786/H/? against metallicity 
for three temperatures: 37 000, 40 000 and 45 000 K. At the 
lower two temperatures, He is partially ionized and thus 
sensitive to changing effective temperature. For the WM- 
basic models, He becomes completely ionized at tempera- 
tures 55 43 000 K. Considering the lower two temperatures 
only, the Kurucz data points increase gradually with abun- 
dance. Since the observed He I A5786/H/3 ratio shows a de- 
crease with increasing abundance (e.g. Bresolin et al. 1999), 
this has led to the suggestion that effective temperature de- 
creases with increasing Z. We see, however, that the WM- 
basic model predictions are not flat but show a decrease 
between 1 and 2 Zq . This is because our new line-blanketed 
models have a softer ionizing flux in the He I continuum 
at high Z (Fig. p"o| ) because of the deeper line-core depths 
in non-LTE. The strength of He I A5876 for a given effec- 
tive temperature therefore decreases with increasing metal 
abundance in a similar fashion to the observations. 

The second significant difference between the WM-basic 
and the Kurucz and CoStar models is the large offset in the 
effective temperature for a given He I line strength at a single 
metallicity. For example at Zq , at a temperature of 37 000 K, 
Hel A5786/H/? = 0.10 (WM-basic); 0.12 (CoStar); and 0.05 
(Kurucz). The absolute effective temperatures therefore de- 
pend heavily on the model atmospheres employed. 



6.2 Cluster models 

Bresolin et al. (1999) compared observational diagnostic di- 
agrams for extragalactic Hn regions with the predictions 
of the cluster models of Leitherer & Heckman (1995). They 
found that the observed line ratios could only be reproduced 
for ages of 1 and 2 Myr, and concluded that either the W- 
R ionizing continua are too hard, or the Hn regions are 
disrupted at a young age. More recently, Bresolin & Ken- 
nicutt (2002) studied a sample of metal-rich extragalactic 
Hn regions. Some of these objects have W-R features in 
their spectra but the nebular line ratios are not affected i.e. 
they find no observational evidence that H II regions ionized 
by clusters containing W-R stars show any signs of harder 
ionizing fluxes. 

In Fig. [if], we plot the "radiation softness" parameter 
r?' defined by Vflchez & Pagel (1988) as 



V = 



[O II] AA3726, 3729/[0 III] AA4959, 5007 



[S II] AA6717, 6731/ [S III] AA9069, 9532 

against the abundance parameter R23 =([0 II] A3727+[0 
III] AA4959, 5007)/H/3. The rf parameter is a sensitive mea- 
sure of the softness of the radiation field and provides an 
excellent test of the correctness of our new models. We show 
the predicted values of rf for metallicities of 2, 1 and 0.2 Z 
and ages of 1-5 Myr (i.e. up to the end of the W-R phase) for 
an instantaneous burst model of mass 10 6 M Q . We have used 



the same photoionization model parameters as in Sect 
and since rf is sensitive to the ionization parameter, we show 
two cases with log (7 = — 3 and —2. We compare the pre- 
dictions of the new grid with the CoStar and SLG92 W-R 
models of SV98 and the observational data of Bresolin et al. 
(1999). It can be seen that the SV98 models are too hard, 
particularly for the log U = —2.0 case, whereas the new 
models are coincident with the data points. We conclude 
that, for the parameter range explored, the ionizing fluxes 
of the new models are in much better agreement with the 
observed emission line ratios of H 11 regions than previous 
model grids. 



7 DISCUSSION 

The success of the technique of evolutionary synthesis for 
deriving the properties of young stellar populations depends 
on the accuracy of the model atmospheres and evolutionary 
tracks for massive stars. In this paper, we have sought to 
provide a grid of realistic ionizing fluxes for O and W-R 
stars using the latest line-blanketed, non-LTE atmosphere 
codes. We have carefully chosen the input stellar param- 
eters to reflect current values available in the literature. 
While these appear to be well established for the O spectral 
class, there have been major revisions for the W-R stars. We 
have incorporated these changes by using lower dumping- 
corrected mass loss rates and lower temperatures compared 
to the previous W-R grid of SLG92. We have also intro- 
duced metallicity-dependent mass loss rates because recent 
observational and theoretical work indicate that the strength 
of W-R winds should scale with heavy metal abundance 
(Crowther et al. 2002a). 

The lower, metallicity-dependent, W-R mass loss rates 
have implications for the amount of mass and energy re- 
turned to the interstellar medium from young starbursts. 
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Figure 15. The predicted strength of the softness parameter rj plotted against the abundance indicator R23 for 2, 1 and 0.2 Zq for 
cluster models with log U = —3 and —2 using (a) the models in this paper and (b) the CoStar and SLG92 models of SV98. The observed 
data points for the extragalactic Hll region sample of Bresolin et al. (1999) are shown for comparison. 



We will address this topic in a future paper. Briefly, in com- 
parison to Leitherer et al. (1992), we find that W-R stars no 
longer dominate the energy input, particularly at low metal- 
licity. Our revised O star mass loss rates, however, are larger 
than those given by the analytical formula of Leitherer et 
al. (1992), and this leads to an overall increase in the rates 
of mass and energy input. 

As discussed in Section ^ and SLG92, the emergent flux 
below 228 A depends on the wind density and hence the scal- 
ing law used to scale the wind parameters with metallicity. 
The precise scaling even for O stars has not been determined 
empirically and relies on the predictions of radiatively-driven 
wind theory. It is of course unknown for W-R stars and we 
simply chose to use the O star scaling law. The predicted 
number of He + ionizing photons at metallicities below solar 
should therefore be considered as rather uncertain. 

We now discuss the presence of nebular Hen A4686 
emission. Guseva et al. (2000) compile a list of 30 Hll 
galaxies with nebular Hen detections. Schaerer (1996) used 
CoStar and the SLG92 models to synthesize the strength of 
nebular Hen A4686 emission in the early stages of a star- 
burst. Schaerer predicts that strong nebular Hell emission 
will be present during the W-R phase with J(Hell)//(H/9) = 
0.01-0.025 for 0.2 Z s: Z sC Z Q . Guseva et al. (2000) find 
that the predicted He n emission line strengths (as given by 
SV98) agree reasonably well with the observations except 
that the SV98 models predict the highest J(Heii)/7(H/3) 
ratios to occur at 1 and 2Zq, while no nebular Hell is de- 
tected in objects with Z > 0.2 Zq. They also find that W-R 
stars cannot be the sole source of nebular He 11 emission since 
only 60 per cent of their sample have W-R stars present, and 



they suggest additional mechanisms such as radiative shocks 
are needed to explain all the observations. 

Assuming case B recombination and an ionization- 
bounded nebula, the strength of He n A4686 relative to H/3 is 
given by J(Heil)//(H/3)= 2.14Q 2 /Q (SV98). If we further 
assume that the minimum detectable 7(Hell)/7(H/3)= 0.01, 
then log Q2/Q0 ^ —2.33 for Hen A4686 to be observable. 
From Fig. hcj, the highest value of logQ2/Qo = —2.4, cor- 
responding to J(Hell)//(H/3)= 0.08, occurs at the begin- 
ning of the W-R phase at 0.2 Zq. Thus any nebular Hen 
will be marginally observable. With the current grid of line- 
blanketed W-R models, it seems unlikely that the observed 
nebular He 11 associated with W-R galaxies is in general due 
to W-R star photoionization. 

We must emphasize that the W-R grid is based on the 
observational parameters of "normal or average" W-R stars 
and thus excludes the rare, hot early WN and WO stars 
with weak winds and significant fluxes above 54 eV. Since 
we do not understand the reasons for their extreme param- 
eters, it is impossible to predict their impact on the ionizing 
output of young starbursts. Given their rarity in the Galaxy 
and LMC, we expect their influence to be negligible. But, 
if for example, they are formed via the binary channel in 
significant numbers at low Z, they may become an impor- 
tant source of He + ionizing photons. We have also neglected 
shocks in the stellar winds; they will increase the ionizing 
flux in the EUV. It is not known how the strength of these 
shocks (caused by instabilities in the line-driving force) will 
scale with metallicity. 

The problem of assigning W-R atmospheres to evolu- 
tionary models is a long-standing one and has been exten- 
sively discussed by SLG92 and in Section 3.2. To connect our 
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new W-R models to the evolutionary tracks, we matched T, 
to 0.6Th y d + 0.4T 2 /3 rather than matching T* directly to 
Tkyd, as in Leitherer et al. (1999) because of the lower tem- 
peratures of the new W-R grid. The differences we find for 
starbursts in the W-R phase compared to the results of SV98 
and Leitherer et al. (1999) is therefore potentially a combi- 
nation of the effects of line blanketing and using lower tem- 
perature W-R models. The overall effect of the method we 
have used to connect the new W-R grid with the evolution- 
ary tracks can be gauged from Figs. ^ and [loj The relative 
heights of the WN and WC bumps at ~ 3.5 and 5 Myr are 
very similar for all three comparisons. This suggests that the 
main differences we find in the W-R phase are primarily due 
to the inclusion of line-blanketing and that the lower tem- 
peratures are a secondary effect. The formalism that we used 
to interface the W-R grid to the evolutionary tracks should, 
however, be considered as rather uncertain. For this reason, 
it is treated as a free parameter in the up-dated Starburst99 
code. 

Finally, we mention that the predictions we have made 
for the W-R phase as a function of metallicity depend on 
the evolutionary tracks of Meynet et al. (1994) for single 
stars with enhanced mass loss rates. Leitherer (1999) has 
compared various W-R synthesis models and discusses the 
differences in single and binary models. He finds that the 
W-R/O and WC/WN ratios are highest when the enhanced 
mass loss tracks are used, and that the W-R phase extends 
to between 6 and 10 Myr when binary formation is included 
(see also SV98). Recent evolutionary models including rota- 
tion (Maeder & Meynet 2000) have removed the necessity 
of artificially increasing the mass loss rate to obtain bet- 
ter agreement with observations of W-R stars. The lifetime 
of the W-R phase is increased, the WC/WN ratio is re- 
duced, and more W-R stars are produced at low Z (Maeder 
& Meynet 2001). The next stage in improving evolution- 
ary synthesis models should be the inclusion of evolutionary 
models with rotation. 



8 CONCLUSIONS 

We have presented a large grid of non-LTE, line-blanketed 
models for O and W-R stars covering metallicities of 0.05, 
0.2, 0.4, 1 and 2Z©. The grid is designed to be used with 
the evolutionary synthesis code Starburst99 (Leitherer et 
al. 1999) and in the analysis of H II regions ionized by single 
stars. We have computed 110 models for O and early B stars 
using the WM-basic code of Pauldrach et al. (2001). The OB 
stellar parameters are defined at solar metallicity and are 
based on the most recent compilations of observational data. 
The mass loss rates and terminal velocities have been scaled 
according to metallicity by adopting power law exponents 
of 0.8 and 0.13 (Leitherer et al. 1992). 

For the W-R grid, we used the model atmosphere code 
CMFGEN of Hillier & Miller (1998) and computed 60 WN and 
60 WC models. The new grid is based on what we consider to 
be the most realistic parameters derived from recent obser- 
vations and individual model atmosphere analyses. In par- 
ticular, the upper temperature limit of the grid is 140 000 K 
since the vast majority of W-R stars that have been anal- 
ysed fall below this value. For the W-R mass loss rates at 
solar metallicity, we used the relationships derived by Nugis 



& Lamers (2000) which are corrected for inhomogeneities 
in the W-R winds. For the first time, we have introduced 
a W-R wind/metallicity dependence and adopted the same 
power l aw e xponents used for the O star models. We argue 
in Sect. |2.2| that recent theoretical and observational work 
indicate that the strengths of W-R winds must depend on 
metallicity. We stress that the new W-R grid excludes the 
few known examples of individual hot W-R stars with weak 
winds that have significant ionizing fluxes above 54 eV. Since 
they are so rare (in the Galaxy and LMC, at least), we ex- 
pect them to make a negligible contribution to the ionizing 
flux of a young starburst. 

We find significant differences in the emergent fluxes 
from the WM-basic models compared to the CoStar mod- 
els of Schaerer & de Koter (1997). For supergiants, the wind 
density determines the transparency below 228 A. Generally, 
we find a lower flux in the He I continuum with important 
implications for nebular line diagnostic ratios. We believe 
that the CoStar models over-predict the number of He ion- 
izing photons through the neglect of photon absorption in 
line transitions and their re-emission at longer wavelengths 
(Crowther et al. 1999). 

We compared the new W-R model emergent fluxes with 
the pure helium W-R models of Schmutz et al. (1992) which 
have very different parameters (particularly wind densities 
and temperatures). We find that at ~ 45 000 K for solar 
metallicity WN models, the emergent flux below 504 A is 
much lower than the SLG92 models because of the inclusion 
of line blanketing. At 60 000 K, we find that blanketing is 
less important, and at 90 000 K, the wind density controls 
the emergent flux below 228 A. 

The WM-basic model grid has been integrated into 
Starburst99 (Leitherer et al. 1999) by directly replacing the 
Lejeune et al. (1997) LTE model library. For the W-R grid, 
we have used a new method of matching the models to the 
evolutionary tracks by taking a weighted mean of the un- 
corrected hydrostatic temperature Tkyd and the corrected 
hydrostatic temperature T2/3. This is necessary because of 
our lower, more realistic temperatures, particularly for the 
WC stars. 

We next compared the output ionizing fluxes of the new 
grid integrated into Starburst99 with the evolutionary syn- 
thesis models of Leitherer et al. (1999) and Schaerer & Vacca 
(1998) for an instantaneous burst and continuous star forma- 
tion. The changes in the ionizing outputs are dramatic, par- 
ticularly during the W-R phase, with the details depending 
on metallicity because of line-blanketing and wind density 
effects. For an instantaneous burst, we find that the number 
of He + ionizing photons (Q2) emitted during the W-R phase 
is negligible at Zq and higher. At lower metallicities, Q2/Q0 
is softer by a factor of ~ 20 compared to the SLG92 models. 
In contrast to Schaerer (1996), we predict that nebular Hell 
will be at, or just below, the detection limit in low metallic- 
ity starbursts during the W-R phase. We also find lower He 
ionizing fluxes for Z ^ 0.4 Zq and ages of 7 Myr com- 
pared to SV98 because of the smaller contributions of the 
line-blanketed models. The ionizing fluxes of the continuous 
star formation models emphasize the differences found for 
the single star models. 

We have tested the correctness of the new model grid 
by computing nebular line diagnostic ratios using CLOUDY 
and single star and evolutionary synthesis models as inputs. 
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For the former, we calculated the nebular He I A5786/H/3 
ratio as a function of Z since this is a sensitive diagnostic 
of Qi/Qo for the temperature range where He is partially 
ionized. Observations indicate that this ratio decreases with 
increasing metal abundance, leading to the suggestion that 
the effective temperature decreases with increasing Z, and 
thus that the upper mass limit may be Z-dependent (e.g. 
Bresolin et al. 1999). The WM-basic models for dwarf O 
stars show a decrease in He I A5786/H/3 above Zq in contrast 
to Kurucz LTE models which are essentially independent of 
Z. This decrease is in the same sense as the observations and 
suggests that the observed decline in He I A5786/H/3 with 
increasing Z is simply due to the effect of line blanketing 
above Zq, and is not caused by a lowering of the upper 
mass limit. 

The ionizing fluxes of the instantaneous burst models 
have been tested by comparing them to the predictions of 
SV98 and the observations of Bresolin et al. (1999) for ex- 
tragalactic Hn regions. In plots of the softness parameter 
rj' against the abundance indicator R23 for ages of 1-5 Myr 
and log U = —2 and —3, we find that the new models cover 
the same parameter space as the data points in contrast to 
the SV98 ionizing fluxes which are generally too hard at all 
ages. 

We therefore conclude that the ionizing fluxes of the 
new model grid of O and W-R stars represent a consider- 
able improvement over the model atmospheres that are cur- 
rently available in the literature for massive stars. To prove 
their worth, they should be tested rigorously against specific 
observations of single Hn regions and young bursts of star 
formation. 

The grid of ionizing fluxes for O and W-R stars 
and the updated Starburst99 code can be obtained from: 
http : //www . star . ucl . ac . uk/starburst . 
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